0D-R148  838  RCTR  RERONRUTICR  ET  RSTRONHUTICfl  SINICfl  (SELECTED 

RRTICLESKU)  FOREIGN  TECHNOLOGV  DIV  WRIGHT-PRTTERSON 
RFB  OH  C  LIU  ET  RL.  27  NOV  84  FTD-ID(RS>T-U58-84 
UNCLASSIFIED  F7G  2874 


TO  FILE  CORY  AD-A148  830 


m 


FTD- ID ( RS )T- 1 1 50-84 


FOREIGN  TECHNOLOGY  DIVISION 


ACTA  AERONAUTICA  ET  ASTRONAUTICA  SINICA 
(Selected  Articles) 


•  -3  -f-'  **  > 


D 'i  iO 

JArJLECTE 


JAN  2  t335 p  §? 
A 


Approved  for  public  release; 
distribution  unlimited. 


84  12  26 


FTD  ~ ID (RS)T-l 150-84 


EDITED  TRANSLATION 


FTD-ID(RS)T-1 150-84  27  November  1984 

MICROFICHE  NR:  FTD-84-C-001116 

ACTA  AERONAUTICA  ET  ASTRONAUTICA  SINICA 
(Selected  Articles) 

English  pages:  38 

Source:  Hangkong  Xuebao,  Vol.  5,  Nr.  1,  1984,  pp.  11-29 

Country  of  origin:  China 
Translated  by:  SCITRAN 

F33657-81-D-0263 
Requester:  FTD/TQTA 

Approved  for  public  release;  distribution  unlimited. 


THIS  TRANSLATION  IS  A  RENDITION  OF  THE  ORIGl- 
NAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 
EDITORIAL  COMNENT.  STATEMENTS  OR  THEORIES 
ADVOCATEDOR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 
ANOOO  NOT  NECESSARILY  REFLECT  THE  POSITION 
OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  Dl- 
VISION. 


PREPARED  BY: 

TRANSLATION  DIVISION 
FOREIGN  TECHNOLOGY  DIVISION 
WP.AFB,  OHIO. 


FTD-ID(RS)T-1150-84 


Date  27  Nov  ]Q  84 


Table  of  Contents 


Graphics  Disclaimer  . . . 

An  Approximated  Solution  of  the  Aircraft  Lateral-Directional  Limit  Cycle 
Oscillation  Induced  by  Aerodynamic  Hysteresis,  by  Liu  Chang  . 

Selection  of  the  Longitudinal  Feedback  Coefficients- of  the  Stability 
Augmentation  System  in  Order  to  Improve  Riding  Qualities,  by  Gao  Hao 

The  Compilation  and  Application  of  a  State-Time  Spectrum  of  Aircraft 
Ambient  Vibration,  by  Gong  Qing  Xiang  . 


Accession  :  £>: 

ntts  grmx 

nic  TAB 


An  Approximated  Solution  of  the  Aircraft  / 

Lateral-Directional  Limit  Cycle  Oscillation  Induced 
by  Aerodynamic  Hysteresis** 

Nanjing  Aeronautical  Institute  Liu  Chang 

ABSTRACT 

This  paper  gives  the  approximated  analytical  equation  and 
determines  the  characteristics  of  the  aircraft  lateral- 
directional  limit  cycle  oscillations  which  are  induced  by 
aerodynamic  hysteresis.  Our  computation  indicates  that  it  agrees 
fairly  well  with  the  result  obtained  by  the  numerical  iteration 
method  in  references  [1,2].  Furthermore,  our  method  eliminates 
the  shortcoming  in  the  references  mentioned  above  in  that  the 
amplitude  of  cycle  oscillation  associated  with  the  roll  rate 
limit  is  not  the  maximum  amplitude. 

I .  INTRODUCTION 

References  [1,2]  pointed  out  that  during  wind  tunnel 
experiments  both  the  roll  in  the  lateral  direction  and  the  yawing 
moment  vary  with  the  yaw  angle  under  the  conditions  of  separation 
of  gas  flow  from  the  surface  of  the  wing  of  the  aircraft  at  high 
angles  of  attack.  Under  these  conditions  there  also  exists  the 
phenomenon  of  hysteresis.  This  model  is  shown  in  Figure  1. 


Figure  1.  Aerodynamic  Rolling  (Yawing)  Moment  Hysteresis 
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The  moment  coefficient  induced  by  yawing  is 

Am, sign(3).  i=x,y  (A) 

a 

In  this  equation  m^S  is  the  moment  coefficient  induced  by  yawing 
when  the  phenomenon  of  hysteresis  does  not  exist,  aok  is  the 
increment  of  the  rolling  moment  coefficient  produced  by 
hysteresis.  This  increment  is  close  to  a  sudden  instantaneous 
change,  comparable  to  a  relay  effect.  Its  magnitude  is  a 
constant  and  its  direction  varies  with  the  sideslip  angle 
velocity  3.  By  using  the  state  variable  analysis  method  of  the 
control  theory  and  under  the  condition  of  neglecting  the  side 
force  equation  from  the  equation  of  motion  in  the  direction 
lateral  to  the  motion,  references  [1,2]  obtained  the  period  and 
amplitude  of  the  lateral-directional  limit  cycle  oscillation 
based  on  a  numerical  iteration  method.  By  using  the  linealized 
harmonic  wave  method  of  the  control  theory,  this  paper  linealized 
the  nonlinear  relay  type  moment  increment  produced  by  hysteresis. 
This  paper  then  applied  the  equilibrium  method  of  harmonic  waves 
[3]  to  obtain  the  analytical  representative  equations  of  the 
frequency  and  amplitude  of  the  limit  cycle  oscillation  based  on 
the  simplified  lateral  equation  of  motion.  From  the  approximated 
equation,  the  effects  of  the  moment  increment  induced  by 
hysteresis  and  the  related  aerodynamic  derivatives  are  very 
obvious.  This  is  very  convenient  to  the  analysis  of  the 
characteristics  of  the  limit  cycle  oscillation.  The  result  of 
the  computation  indicates  that  the  result  given  by  the 
approximated  analytical  equation  is  in  good  agreement  with  the 
result  given  by  references  [1,2]. 

The  rules  of  the  aerodynamic  manual  [4]  will  be  followed  for 
all  the  symbols  which  are  used  in  this  paper  but  not  explained 
explicitly.  Further  discussion  on  the  related  references  [1,2] 
can  be  found  in  reference  [5]. 
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II.  COMPUTATIONAL  METHOD 

We  will  consider  the  side  force  equation  for  the  lateral 
perturbating  notion  of  the  aircraft.  We  will  also  assume  that 
the  attack  angle  and  all  control  surfaces  remain  unchanged. 


In  this  way  we  can  obtain  im  dt  mM*&+  \f~'u. 


If  the  aircraft  maintains  a  straight  flight  path, 
df/dt=d&/dt=6,  and  the  equation  shown  above  can  be  simplified 
into  ,  , 


+  WT»3  -r 


The  derivatives  in  this  equation  have  short  lines  on  top  of  them, 

Indicating  that  these  derivatives  are  divided  by  their 

corresponding  rolling  moments,  such  as  My/Iy,  Mx=Mx^x 

The  nonlinear  additional  moment  coefficient  ’m  is  induced 

by  aerodynamic  hysteresis  is  the  relay  type,  as  shown  in  Figure 

2.  We  can  obtain  the  equivalent  moment  coefficient  derivative  by 

[31 

using  the  linearized  harmonic  wave  method  .  It  can  be 
represented  as : 


Or  we  can  also  represent  it  as  the  moment  derivative 


In  this  equation  aM^  *  AM^qsl/I^,  i  =  x,y  and  0^,  u»^  are  the 
amplitude  and  frequency  of  the  limit  cycle  oscillation. 

By  substituting  the  increment  of  the  hysteresis  moment 
derivative  (2)  into  equation  (1),  we  can  obtain  the  equation 
which  considers  the  equation  of  lateral  perturbating  motion  under 
the  influence  of  mement  hysteresis. 


Figure  2.  Nonlinear  Additional  Moment  Coefficient. 

1.  Rolling  Moment  Hysteresis 

After  taking  the  additional  increment  of  the  rolling 
hysteresis  moment  aM®  into  consideration,  equation  (1)  can  be 
written  in  the  following  form: 


<-st 


(3) 


By  eliminating  a>x  from  equation  (3),  we  can  get 


(4) 


Based  on  the  harmonic  wave  equilibrium  method  [3],  we  can 
assume  the  approximated  solution  of  the  limit  cycle  oscillation 
to  have  the  following  form.  We  then  have 


0  -e,ccsc 0,1 

a—— 


(5) 
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Substituting  equation  (5)  into  equation  (4)  and  after 
rearranging,  we  have 


/ 13 
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By  comparing  the  corresponding  first  order  harmonic  wave 
coefficients  on  two  sides  of  the  equation  above,  we  have 


(6) 


•*®i  —  [mJ+mJ'  («+-££&•) 

0  m‘-(Mi-hM7')9!toi  +  (M7’ia!-M:rM!)6l  .  (7) 


From  equation  (7)  we  can  obtain  the  frequency  of  limit  cycle 
oscillation  at  sideslip  angle  8  as 


<o  mT’m}  y 
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(8) 


By  substituting  equation  (8)  into  equation  (6),  we  can 
obtain  the  corresponding  amplitude  of  the  limit  cycle  oscillation 
as  _ _ 

(9) 


SOI ,  (to J + «}  +  J  -  mVmI  )  » 

From  the  above  equation,  it  is  clear  that  the  magnitude  of 
the  amplitude  of  limit  cycle  oscillation  is  directly  proportional 
to  the  additional  rolling  moment  AMx  induced  by  the  aerodynamic 
hysteresis. 

From  the  first  equation  in  equation  C 3 ) ,  we  can  obtain 


®.--=rr  (8-m?8-mJ3). 

M  i 


(10) 


We  now  substitute  equation  (5)  into  equation  (10),  after 
rearranging  the  equation  we  can  get 


Cj5f;9,ci),sjn«,<— (0i»?  +  m59 Jcofiu;0.  (E) 
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the  above  equation 


Let  C  =  ,  C^=  -(^  + 

can  be  rewritten  as: 


C,  -  m  Je ;  <» ,  /mV,  C.  -  -  ( 8 .  o>|  +  m  i  8 ; )  /  M?'. 

'flD ,  4«  C ,  si  n  a ,/  +  C  .cos  ca . =*  0 ;  s :  n  ( <o  ,<  +  <? , ) 

In  this  equation  we  have  __  __ 

0.-i/  CJ  +  CJ 
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(11) 


(F) 


According  to  equation  (11)  we  can  determine  the  amplitude  of 
the  limit  cycle  oscillation  e2  at  the  rolling  speed  as 


9. 


(mt ‘»Dl+(w{+A(;)i . 


(12) 


From  the  equation  above  we  can  see  that  the  oscillating 
amplitude  e2  at  rolling  speed  u»x  is  linearly  proportional  to  the 
oscillation  ex  at  sideslip  angle  8. 

2.  Yawing  moment  hysteresis 

After  taking  into  consideration  the  additional  increment  of 
the  yawing  hysteresis  moment  derivative  aM£ ,  equation  (1)  can  be 
rewritten  into  the  following  form 


(13) 


(14) 


By  eliminating  lii^from  equation  (13),  we  can  get 
3~(  +  -m?* 

Based  on  tne  same  principle  and  according  to  the  harmonic 
wave  equilibrium  method,  the  frequency  u>2  of  the  limit  cycle 
oscillation  at  sideslip  angle  S  can  be  obtained  from  the  equation 
above 


At  =*  -~A.W, 

-  \tV\fr  (G) 


and  the  amplitude  of  the  limit  cycle  oscillation  6, 


_ 4  _ 

J(/4,coj-C,) 


(16) 


We  can  see  from  equation  (16)  that  the  magnitude  of  the 
amplitude  9,  is  directly  proportional  to  the  increment  of  the 
yawing  hy teres is  moment. 

Similarly,  after  taking  the  additional  increment  of  the 

•a 

yawing  hysteresis  moment  derivative  into  consideration,  we  can 
determine  the  amplitude  04of  the  limit  cycle  oscillation  at 
rolling  rate 


9, 


_9i_ 
siV  V 


(17) 


From  the  equation  above  we  can  see  that  the  amplitude  of 
oscillation  0*  at  rolling  rate  u  is  directly  proportional  to  the 
amplitude  of  oscillation  0,  at  sideslip  angle  S  . 

3.  The  Amplitude  and  Frequency  at  Rolling  Angle 
According  to  the  approximated  relationship  equation  dy/dt=(u 
we  can  determine  the  amplitude  of  the  limit  cycle  oscillation 
at  rolling  angle  T  as 


r 


In  this  equation  and  <u  are  the  amplitude  and  frequency  of  the 
rolling  rate  limit  cycle  oscillation. 
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Ill  .  SAMPLE  CALCULATION  AND  ANALYSIS 

By  taking  the  data  set  corresponding  to  the  landing  state  of 
the  F-94A  aircraft  as  reference  [1,23  did,  we  can  compute  the 
corresponding  moment  derivatives  as: 

M®  =  -2.8256  1/second  =  -1.3214  1/second 

Mux  =  -2.4593  1/second  M*"  _  =  0.06283  1/second 

M^  =  -1.5193  1/second  M^  =  =  -0.2489  1/second 

x  y  y 

According  to  the  equations  derived  earlier,  we  can  compute 
the  amplitude  9  and  period  T  (  =  — ^— )  of  the  rolling  and  yawing 
lateral-directional  cycle  oscillation  induced  by  hysteresis.  The 
results  are  shown  in  Table  1.  For  easy  comparison,  we  have  also 
listed  the  results  obtained  by  the  numerical  iteration  method  of 
reference  [23  in  Table  1. 
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1.  Table  1.  Amplitudes  and  Periods  of  Sideslip  Angles  and  Roll 
Rate  Limit  Cycle  Oscillations 

2.  type  of  hysteresis 

3.  magnitude  of  hysteresis 

4.  9  (radian) 

5.  9b  (1/radian) 

6.  Tvsecond) 

7.  according  to  our  paper's  analytical  equation 

8.  according  to  the  numerical  iteration  method  of 
reference  [23 

9.  according  to  our  paper's  analytical  equation 

10.  according  to  the  numerical  iteration  of 
reference  [23 

11.  according  to  our  paper's  analytical  equation 


12.  according  to  the  numerical  iteration  method  of 
reference  C2] 

13.  rolling 

14.  yawing 

From  this  table  we  can  see  that: 

(1)  The  period  of  the  limit  cycle  oscillation  T  computed 
based  on  the  analytical  equation  proposed  by  this  paper  is  very 
close  to  the  period  obtained  by  using  the  numerical  iteration 
method  of  reference  C2]. 

(2)  With  the  exception  of  the  0.,  induced  by  rolling 
hysteresisAM  ,  the  amplitude  of  the  limit  cycle  oscillation 

X 

computed  based  on  the  analytical  equation  proposed  by  this  paper 
is  also  in  good  agreement  with  the  computed  one  based  on  the 
numerical  iteration  method  of  reference  [2].  Reference  [2]  gave 
the  computer  simulated  result  of  the  lateral  equation  of 
perturbating  motion  (1)  which  takes  into  consideration  the  moment 
increment  induced  by  hysteresis.  The  curves  of  the  transitional 
processes  are  shown  in  Figures  3  and  4.  The  results  show  the 
reliability  of  the  approximation  method. 

(3)  The  amplitude  of  the  rolling  rate  limit  cycle 
induced  by  the  rolling  hysteresis  as  computed  by  the  method  of 
this  paper  (,<*0.42)  is  higher  than  the  value  computed  by  the 
numerical  iteration  method  of  reference  [2]  (0.29).  But  we  can 
see  from  Figure  3  that  it  is  relatively  close  to  the  value 
obtained  by  an  integration  of  the  original  numerical  value  of  the 
equation  (**0.38).  Actually,  the  g^  value  obtained  from  the 
numerical  iteration  method  of  reference  [2]  is  a  value 
corresponding  to  the  condition  where  b  vanishes  to  zero.  This 
corresponds  to  the  lower  peak  value  in  Figure  3,  not  the  maximum 
amplitude  of  the  rolling  rate  limit  cycle  oscillation.  Of 
course,  the  computational  result  of  the  analytical  equation  will 
inevitably  bring  in  a  certain  amount  of  error  since  the  rolling 
rate  reaction  differs  from  the  wave  form  of  the  first  order 
harmonic  wave  significantly  under  certain  conditions. 


Figure  3.  Limit  Cycle  Oscillation  Induced  by  Rolling 

Hysteresis 

1.  a*  (radian /second)  8  (radian) 

2.  9w  computed  according  to  this  paper 

3.  0«x  computed  according  to  reference  [2] 

4.  t(§econd) 

Similarly,  we  have  also  conducted  computations  by  using  the 
e«ox  and  T  values  in  Table  1  in  the  approximation  equation  (18) 
which  determines  the  rolling  angle  r  and  amplitude  0r.  The 
results  obtained  are  also  relatively  satisfactory. 


Figure  4.  Limit  Cycle  Oscillation  Induced  by  Yawing  Hysteresis 

1.  a*  (radian /second)  g  (radian) 

2.  t“( second) 

IV.  CONCLUSION 

The  approximated  analytical  equation  in  this  paper  can  be 
used  to  determine  the  characteristics  of  the  lateral-directional 
limit  cycle  oscillation  induced  by  aerodynamic  hysteresis.  The 
results  of  the  computation  indicate  that  they  are  in  good 
agreement  with  the  results  obtained  by  using  the  numerical 
iteration  method  of  references  Cl, 2].  This  paper  also  avoids  the 
shortcoming  of  these  references  in  that  the  amplitude  of  the  roll 
rate  limit  oscillation  induced  by  rolling  hysteresis  is  not  the 
maximum  amplitude.  In  addition  to  the  simplicity  in  computation, 
the  analytical  equation  is  convenient  for  the  analysis  of  the 
effects  of  the  magnitude  of  hysteresis  and  the  various  related 
aerodynamic  derivatives  on  the  characteristics  of  the  limit  cycle 
oscillation  since  the  equation  is  very  explicit. 

Professor  Tan  Zhen  Hua  has  provided  many  valuable  comments. 
We  wish  to  express  our  appreciation  here. 
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Abstract 

Wind  tunnel  tests  on  aircraft  models  often  demonstrate  that  the  vari¬ 
ation  of  lateral-directional  rolling  and  yawing  moments  with  the  yaw 
angle  gives  evidence  of  aerodynamic  hysteresis  occurred  at  high  angle  of 
attack- The  non-linear  additional  increments  of  the  moment  due  to  aero¬ 
dynamic  hysteresis  are  linearized  by  the  harmonic  linear  method.  Then, 
based  on  the  harmonic  balance  method,  an  approximate  analytical  formula 
is  derived  for  determining  the  characteristics  of  aircraft  lateral-directional 
limit  cycle  oscillation  induced  by  aerodynamic  hysteresis.  Calculations  by 
this  method  show  that  the  results  are  in  good  agreement  with  data  ob¬ 
tained  by  the  numerical  iterating  method  in  referencesCl)andC2).  The  me¬ 
thod  eliminates  the  shortcoming  in  the  above  references  that  the  ampli- 
,  tude  of  roll  rate  limit  cycle  oscillation  induced  by  roil  hysteresis  is  not  the 
maximum. 


Selection  of  the  Longitudinal  Feedback  Coefficients 
of  the  Stability  Augmentation  System  in  Order  to 
Improve  Riding  Qualities** 

Northwestern  Polytechnical  University  Gao  Hao 

ABSTRACT 

The  longitudinal  feedback  coefficients  of  a  stability 
augmentation  system  (SAS)  for  a  sample  aircraft  were  selected 
based  on  the  condition  that  the  normal  acceleration  response  at 
the  pilot's  position  was  approximately  zero.  The  improvement  in 
riding  qualities  was  estimated  based  on  references  [2.3].  In 
addition,  the  dynamic  characteristics  of  the  combined  aircraft- 
SAS  system  were  also  analyzed.  The  results  indicate  that  the 
requirement  of  riding  quality  could  be  satisfied  when  the  control 
anticipative  parameter  (CAP)  of  the  combined  aircraf t-SAS  system 
was  close  to  9.81 /lx,  and  that  the  working  range  of  the  SAS  was 
within  10%  of  the  maximum  inclination  angle  of  the  flat  tail.  If 
the  area  of  the  flat  tail  was  relatively  large  compared  to  the 
area  of  the  wing,  such  as  when  Ybz/7°*  0.2,  the  lift  term  due  to 
the  elevator  deflection  would  have  to  be  considered  in  the 
analysis  of  the  dynamic  characteristics  of  the  combined  aircraft- 
SAS  system. 

INTRODUCTION 

For  aircrafts  which  will  have  to  take  low  altitude  combat 
missions,  the  design  of  the  exterior  arrangement  of  the  aircraft 
alone  is  generally  not  going  to  be  enough  to  satisfy  the 
requirement  of  riding  quality.  The  results  of  our  analysis  and 
computation  indicate  that  we  should  try  to  improve  the  riding 
quality  of  the  aircraft  mainly  by  trying  to  reduce  the  normal 
acceleration  response  at  the  pilot's  position  when  the  aircraft 
is  flying  in  atmospheric  turbulence.  In  order  to  achieve  this, 
we  can  introduce  various  types  of  automatic  devices  in  the 
longitudinal  passage  of  the  aircraft.  The  effects  of  the  various 
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types  of  automatic  devices  on  the  improvement  of  riding  quality 
can  be  found  in  reference  E 3 ] -  By  using  a  certain  type  of 
intercepting  fighter  as  an  example,  this  paper  analyzed  several 
problems  involving  the  selection  of  the  feedback  coefficients  of 
the  stability  augmentation  system  (SAS)  using  the  improvement  in 
riding  quality  of  the  aircraft  by  the  longitudinal  SAS  as  a 
guide.  The  standards  of  the  aircraft  riding  qualities  used  in 
this  paper  followed  the  related  specification  of  the  United 
States  Military  manual  MIL-F-8785B. 

I.  The  Selection  of  the  Feedback  Coefficients  of  the 
Longitudinal  Augmentation  System  Used  to  Improve  Riding  Quality 

In  principle,  we  can  satisfy  the  standard  of  riding  quality 
by  using  various  ways  to  select  the  feedback  coefficients  of  the 
stability  augmentation  system.  Here  we  will  only  discuss  the 
selection  of  the  feedback  coefficients  of  centered  SAS  based  on 
the  criterion  that  the  requirement  of  riding  quality  be 
satisfied. 

If  we  can  use  the  short  period  equation  of  motion  to 
approximately  describe  the  dynamic  characteristics  of  the 
aircraft,  and  we  can  take  the  diving  angular  velocity  of  the 
aircraft  reaction,  the  attack  angle,  and  the  increment  of  the 
loading  coefficient  in  the  normal  direction  as  the  feedback 
signals,  the  regulatory  principle  of  the  SAS  under  ideal 
conditions  is 
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The  block  diagram  of  the  combined  aircraft-SAS  system  is  shown  in 
Figure  1. 

In  order  to  simplify  the  discussion,  we  will  not  consider 
the  dynamic  characteristics  of  the  auxiliary  propeller  and 
steering  mechanism.  This  is  to  say  that  we  will  set  Gi*G2“  1. 
From  this  we  can  transform  Figure  1  into 
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Figure  1.  Block  diagram  of  the  aircraft-SAS  System. 

In  this  figure,  G^,  G£  are  the  transmission  function  of 
the  auxiliary  propeller  and  steering  mechanism, 
respectively; 

Enlargement  coefficient  and  transmission  function 
which  are  related  to  the  force  of  the  pilot  control  stick; 
Ka,  K*.  Kn  -  the  feedback  coefficient  of  the  attack  angle, 
the  inclination  angle,  and  the  loading  coefficient; 
m,  /  ft.  *  Ad'/Adu  -transmission  function  of  the  aircraft. 
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In  this  equation  we  have 
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From  this  we  can  see  that  the  introduction  of  the  fl*  A*y,Aa 
signals  is  equivalent  to  a  change  in  the  frequency  of  the  short 
period  mode  of  the  aircraft,  while  only  the  introduction  of  the  S 
signal  can  change  the  damping  of  the  short  periodic  state. 

In  principle,  we  can  select  any  arbitrary  combination  of  K  , 
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K  and  K  to  satisfy  the  requirement  of  the  short  period  mode 
*  n 

described  above.  Here  we  will  start  our  analyses  from  the 
requirement  of  the  riding  quality  of  the  aircraft.  This  is  to 
say  that  we  will  start  from  the  condition  that  the  normal- 
directional  loading  coefficient  reaction  at  the  pilot's  position 
be  kept  at  a  minimum. 


We  will  assume  that  under  the  action  of  exterior 
perturbation  the  increment  of  the  normal-directional  loading 
coefficient  reaction  at  the  center  of  mass  of  the  aircraft  can  be 
represented  by  an^,  the  angular  acceleration  can  be  represented 
as  »  ,  and  the  distance  between  the  pilot's  seat  and  the  center 
of  mass  of  the  aircraft  is  lx<  The  condition  which  guarantees 
the  normal-directional  loading  reaction  at  the  pilot's 
position  to  be  zero  is 


After  simplification  and  some  transformation  which  is  not 
complicated,  the  condition  above  can  be  rewritten  as 
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Generally  speaking,  1  are  smaller  than  9.81  meters  for 
fighters,  attack  planes  and  figther  bombers.  This  is  why  we  can 
expect  the  CAP  values  obtained  from  the  equation  above  to  satisfy 
the  standard  of  riding  quality. 

By  expanding  the  u  and  (x^/a),  s-«  terms,  we  can  ultimately 
obtain  the  condition  which  determines  CAP  as 
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For  simplicity,  we  can  set  K*  and  to  have  a  constant 
proportionality  relationship.  In  other  words 
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The  value  of  K  can  be  determined  through  a  sample  calculation. 
At  the  same  time,  we  believe  that  the  introduction  of  the  K 
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signal  primarily  compensates  for  the  effect  of  the  variation  in 
the  intrinsic  frequencies  of  the  short  period  modes  of  the 


aircraft  when  the  aircraft  is  flying  under  different  M  numbers. 
Because  of  this,  we  can  select  KaQ  by  using 

(H) 

Ultimately  we  can  obtain  the  equations  for  the  selection  of  K„ 

and  K  as  follows: 
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Reference  [4]  recommends  0.02<  K<  0.08.  The  results  of  the 
sample  calculation  indicate  that  the  riding  quality  is  better 
when  K  is  close  to  the  higher  limit.  In  addition,  according  to 
the  equation  above  the  three  feedback  coefficients  should  be 
functions  of  the  flight  altitude  and  the  M  number.  The 
calculation  indicates,  however,  that  they  do  not  vary  a  lot 
within  a  flight  range  of  M<1  so  they  can  be  set  as  constants. 

II.  SAMPLE  CALCULATION 

In  order  to  verify  that  the  feedback  coefficients  we  have 
selected  were  reasonable,  we  have  conducted  a  sample  calculation 
by  using  a  certain  type  of  fighter  plane  as  an  example.  We  will 
now  explain  some  of  the  results  of  the  calculation. 

1.  Figure  2  indicates  the  variations  of  the  feedback 
coefficients  with  the  flight  altitude  and  the  M  number.  It  can 
be  seen  from  the  figure  that  within  the  range  of  M<1.0,  the  K  , 
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and  varied  only  slightly  with  the  altitude  and  the  M  number 
for  the  aircraft  used  in  the  sample  calculation.  We  can  set  them 
to  be  constants  in  practical  applications. 
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1.  aH  a  0.05  Km 

2.  H*  1 IKm 
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Figure  2.  Variation  of  the  Feedback  coefficients  with 
the  Mach  number 

2.  After  the  three  feedback  coefficients  are  set  to  be 

constants,  an  example  of  the  variation  in  the  characteristics  of 

the  short  period  mode  corresponding  to  the  combined  aircraft-SAS 

system  is  shown  in  Table  1.  It  can  be  seen  from  this  table  that 

a  reasonable  selection  of  the  K  ,  Kt,  and  K  can  increase  the 
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relative  damping  ratio  and  the  intrinsic  frequency  at  short 
periods,  hence  increasing  the  stability  at  short  periods.  At  the 
same  time,  the  sensitivity  of  the  aircraft  to  the  control 
reaction  increases  due  to  the  increase  in  the  intrinsic  frequency 
and  the  CAP  value.  As  for  the  effects  of  the  feedback 
coefficients  on  the  other  control  characteristics,  they  can  be 
obtained  by  the  reasonable  selection  of  Kf  and  G.  We  will  not 
discuss  this  aspect  of  the  problem  in  the  present  paper. 

3.  The  effect  on  the  long  period  state  by  the  introduction 
of  the  stability  augmentation  system.  Table  2  lists  the 
variations  of  the  various  coefficients  in  the  characteristic 
equation  before  and  after  the  combination  of  the  sample  aircraft 


and  the  SAS.  It  can  be  seen  that  under  the  two  kinds  of  flight 
conditions  in  Table  2,  a^<0  when  the  SAS  is  not  installed  and 
the  long  period  mode  is  unstable.  After  the  SAS  is  installed, 
a^  >0  and  R>0  and  the  long  period  mode  is  stable. 

4.  The  effect  of  the  SAS  on  the  riding  quality  of  the 
aircraft.  We  have  computed  the  following  indices  according  to 
references  C  2 , 3  3 : 
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1.  Table  1.  The  effect  of  the  stability  augmentation 

system  on  the  characteristics  of  the  short 
period  mode 

2.  (kilometer) 

3.  condition 

4.  category 

5.  without  the  stability  augmentation  system 

6.  oihu  (1  /radian) 

7.  with  the  stability  augmentation  system 

8.  w  (1/radian) 
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r  1  Table  2.  The  effect  of  the  stability  augmentation 

■  *  system  on  the  characteristics  of  the  long 

period  mode 
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2.  (H=0 .05  Km) 

3.  without  the  stability  augmentation  system 

4.  with  the  stability  augmentation  system  added 

5.  without  the  stability  augmentation  system 

6.  with  the  stability  augmentation  system  added 

7.  a»(1 /second^ 

8.  aa (1/secondj ) 

9.  a,  ( 1 /second* ) 

10.  a* (1 /second  ) 

(1)  Within  each  minute  of  the  flight  time,  we  will  set  the 
number  of  times  in  which  the  root  means  square  value  of  the 
reaction  of  the  normal  directional  acceleration  of  the  pilot's 
position  exceeds  0.5g  to  be  N.  The  relationship  between  N  and 
the  evaluation  of  the  passenger  are  as  follows: 

N  ^5  comfortable 
5<  N  ^13  acceptable 
N  >13  not  acceptable 

(2)  The  passenger  work  productivity  index  Hg.  Reference  _ 
C23  believes  that  He  should  be  below  0.035.  It  believes  that  Hg 
0.045  is  unacceptable. 
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The  computational  results  of  the  above  mentioned  Indices 
for  the  sample  aircraft  are  listed  in  Table  3. 
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1.  Table  3.  The  Riding  Quality  of  the  Sample  Aircraft 

2.  (Km) 

3.  Stability  augmentation  system 

4.  without 

5.  with 

6.  index 

7.  remark 

8.  should  be  <5 

9.  should  be  <0.045 

10.  should  be  <5 

11.  should  be  <0.045 

In  this  table,  the  following  equation  is  used  to  compare 
the  reaction  of  the  normal-directional  acceleration  at  the 
pilot's  position 
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When  the  input  quantities  are  selected  as  the  vg>ag»  and 
wx  induced  by  the  atmospheric  turbulent  flow,  the  frequency 
splctrum  of  Aan  is 
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In  this  equation  G.^,  3  used  the  Vg ,  ag,<uxg  as  input 

quantity,  respectively.  Aan  is  the  transmission  function 
corresponding  to  the  output  quantity.  The  footnote  *  indicates 
the  complex  conjugate. 

The  numerical  computation  indicates  that  *  is  primarily 
determined  by  the  second  term  in  the  equation  above. 

It  can  be  seen  from  Table  3  that  the  riding  quality  of  the 
original  aircraft  was  relatively  poor  under  the  computed  flight 
state.  The  requirement  of  the  riding  quality  index  was  bascially 
staisfied,  however,  after  the  stability  augmentation  system  was 
introduced.  This  was  due  to  the  fact  that  the  stability 
augmentation  system  increased  the  intrinsic  frequencies  at  the 
center  of  energy  distribution  of  the  atmospheric  turbulence.  At 
the  same  time  the  stability  augmentation  system  also  increased 
the  relative  damping  ratios  at  short  periods.  All  these 
contributed  to  reduce  the  reaction  of  the  normal-directional 
acceleration  when  the  aircraft  was  flying  within  the  atmospheric 
turbulence,  hence  improving  the  riding  quality.  From  these 
points  we  can  see  that  the  method  used  by  this  paper  to  select 
the  feedback  coefficients  of  the  longitudinal  stability 
augmentation  system  can  be  applied  in  principle  to  the  stability 
augmentation  systems  with  the  primary  objectives  of  improving  the 
riding  qualities. 

5.  Reference  [4]  recommends  that  the  working  range  (limit) 
of  the  stability  augmentation  system  should  be  less  than  107.  of 
the  maximum  inclination  angle  of  the  flat  tail.  All  the 
aircrafts  used  in  the  computation  are  within  this  range. 

6.  The  variation  in  the  dynamic  characteristics  of  the 
combined  aircraft-stability  augmentation  system  within  the 
contribution  of  the  lift  induced  by  the  flat  tail  is  considered. 

As  a  result  of  the  design  of  the  exterior  arrangement  of  the 
modern  fighters,  it  is  possible  for  the  area  of  the  flat  tail  to 
be  close  to  207.  of  the  area  of  the  wing.  At  this  time  we  have 
to  consider  the  variation  in  the  lift  induced  by  the 


inclination  of  the  flat  tail,  which  is  the  *Y^ZA6  term.  After 
taking  this  term  into  consideration,  the  short  period 
transmission  function  of  the  aircraft  is 
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By  comparing  these  results  with  the  results  which  did  not  take 

the  flat  tail  lift  into  consideration,  we  found  that  there  is 

only  a  numerical  quantity  difference  for  u>  /A6  .  For  Aa/ a6  , 

there  is  an  additional  zero  point  in  addition  to  a  difference  in 

quantity.  There  is  an  additional  second  order  link  on  the 

numerator  of  Ax  /A6_  and  the  characteristic  equation  remains 

y  z 

unchanged.  The  computational  results  of  the  sample  calculation 
indicate-,  that  when  Y^z  /  Ya  0.2 Q  the  above  mentioned 

variation  will  induce  relatively  large  variations  in  the  dynamic 
characteristics  of  the  combined  aircraft-SAS  system  so  it  cannot 
be  ignored. 

III.  CONCLUSIONS 

1.  For  the  longitudinal  stability  augmentation  system  with 
the  main  objective  of  improving  the  riding  quality  of  the 
aircraft,  we  can  select  the  feedback  coefficients  of  the 
stability  augmentation  system  by  using  the  approximated  condition 
that  the  reaction  of  normal-directional  acceleration  at  the 
pilot's  position  be  kept  at  a  minimum.  This  type  of  combined 
aircraft-SAS  system  can  satisfy  the  requirement  of  flight  quality 
at  both  short  and  long  periods,  and  it  also  assures  that  the 
aircraft  has  a  good  riding  quality. 


2.  Within  a  range  where  the  M  number  is  not  large,  such  as 
M  <1,  the  feedback  coefficients  selected  according  to  the 
procedures  can  be  set  as  constants , and  they  also  do  not  vary  with 
the  flight  altitude. 

3.  If  the  area  of  the  flat  tail  is  large  relative  to  the 

— b  —a 

area  of  the  wing,  such  as  when  Y  z  /y  ^  0.20,  we  then  have  to 

consider  the  Y^z  term  for  the  analysis  of  the  dynamic 
characteristics  of  the  combined  aircraft-SAS  system. 
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SELECTION  OF  SAS'  LONGITUDINAL  FEEDBACK 
COEFFICIENTS  TO  IMPROVE 
RIDING  QUALITIES 

Gao  Hao 

( Northwestern  Polytechnical  University ) 

Abstract 

The  longitudinal  feedback  'coefficients  of  a  stability  augmentation 
system  (SAS)  for  a  sample  aircraft  were  selected  in  the  light  of  the  con¬ 
dition  that  the  normal  acceleration  response  at  the  pilot's  position  appro¬ 
ximated  to  zero.  With  the  aid  of  references  C22  and  C-3  }  the  improve¬ 
ment  in  riding  qualities  of  the  airframe/SAS  system  due  to  the  selected  co¬ 
efficients  were  evaluated.  The  dynamic  characteristics  of  the  airframe/SAS 
system  were  also  analyzed. 

The  results  demonstrate  that  the  riding  qualities  are  desirable  when 
the  control  anticipative  parameter  of  the  airframe/SAS  system  approaches 
to  9.8 l/l„  and  the  operating  range  of  the  SAS  is  limited  to  the  elevator 
displacement  less  than  10%  of  <the  max..displacement.  It  is  necessary  to 
take  into  account  the  lift  due  to  the  eievator  deflection  in  the  analysis 
of  the  dynamic  characteristics  of  the  airframe/SAS  system*  if  the  tail 
area  is  large  enough  in  comparison  with  the  wing  area  such  as 
0.2. 
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The  Compilation  and  Application  of  A  State-Time 
Spectrum  of  Aircraft  Ambient  Vibration* 

Nanchang  Aircraft  Manufacturing  Company  Gong  Qing  Xiang 

ABSTRACT 

The  state-time  spectrum  of  the  aircraft  vibration  is  an 
important  statistical  parameter  of  the  aircraft.  It  is  a 
prerequisite  for  inducing  the  ambient  vibration  data  of  an 
aircraft  with  the  state  inductive  method.  This  paper  proposes  a 
compiling  method  for  the  state-time  spectrum  of  aircraft  ambient 
vibration.  This  method  is  based  on  the  foundation  of  a  great 
volume  of  flight  statistical  data,  and  it  also  fully  considers 
both  the  flight  condition  in  the  air  as  well  as  the  condition  of 
the  ground  tests.  The  state-time  spectrum  of  aircraft  ambient 
vibration  compiled  this  way  can  better  represent  the  condition  of 
the  practical  application  of  the  aircraft.  The  compiling  method 
has  already  been  used  in  engineering  applications. 

SYMBOLS 

A  the  average  ratio  between  the  flight  time  and  the  time 
of  ground  tests  per  each  flight 
B  the  number  of  ground  tests  within  the  life  span  of 
the  aircraft 

c  the  number  indicating  the  state  of  the  ground  test 

m  the  number  indicating  the  flight  course 

1  number  of  aircrafts  being  included  in  the  statistics 
the  actual  take  off  landing  numbers  of  the  i-th 
aircraft  for  the  j-th  course 
Nj  the  total  actual  take  off  landing  numbers  of  all 
the  aircrafts  included  in  the  statistics  for  the 
j-th  course 

N  the  total  actual  take  off  landing  numbers  of  all 
the  aircrafts  included  in  the  statistics 
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the  average  take  off  landing  numbers  during  the 
application  life  span  of  an  aircraft 
the  take  off  landing  numbers  for  a  K-type 
typical  course  during  the  application  life  span 
of  an  aircraft 

the  take  off  landing  numbers  of  an  aircraft 
within  the  application  life  span  of  the 
equipment  loaded  in  the  aircraft 
number  of  typical  courses 
number  of  typical  states 

the  flight  time  of  the  i-fch  aircraft  for  the 
j-th  course 

the  total  flight  time  of  all  the  aircrafts 
included  in  the  statistics  for  the  j-th  course 
the  total  flight  time  of  all  the  aircrafts 
included  in  the  statistics 
the  average  time  duration  per  flight 
the  time  used  for  the  K-th  type  typical  course 
by  all  the  aircrafts  included  in  the  statistics 
the  application  life  span  of  the  aircraft  on 
the  onboard  equipment 

the  vibration  life  span  of  the  onboard  equipment 
the  flight  time  used  for  the  K-th  type  typical 
course  during  the  application  life  span  of  the 
aircraft 


the  time  occupied  by  the  Q-th  type  state  during 
the  flight  for  the  K-th  type  typical  course 
the  time  occupied  by  the  Q-th  type  flight 
state  during  the  application  life  span  of  the 
aircraft 

the  average  time  occupied  by  the  Q-th  type 
flight  state  per  each  flight 


t  the  average  time  taken  by  the  a-th  ground 
01 

test  state  per  each  flight 

the  time  taken  by  the  a-th  ground  test  state 
during  the  application  life  span  of  the 
aircraft 

t^  the  average  time  necessary  to  conducted  an 
extended  vibration  experiment  per  each 
flight 

T  the  time  used  for  the  extended  vibration 
e 

experiments  for  the  onboard  equipment  with 
an  application  life  span  of  T 

s 

t  the  average  time  used  for  each  ground  test 

t^  the  time  used  for  the  k-th  type  typical 

course  per  flight 

Vx  the  average  time  taken  by  the  x-th  ground 

test  state  per  each  ground  test  process 

Wq  the  value  of  random  vibration  corresponding 

to  the  Q-th  type  typical  state 

Vf  „  the  maximum  value  of  random  vibration  of  all 
max 

the  typical  states 
I .  INTRODUCTION 

In  order  to  specify  the  experimental  standards  of  the 
onboard  equipment  and  the  vibration  of  the  instruments  it  is 
necessary  to  conduct  vibration  experiments  on  the  aircraft 
Then  we  can  proceed  to  analyze  the  original  vibration  data 
measured  during  the  experiments.  In  addition,  it  is  necessary 
for  us  to  use  the  appropriate  statistical  methods  to  synthesize 
and  induce  the  large  amount  of  analytical  data.  Only  by  doing 
this  can  we  come  up  with  the  basic  experimental  parameters 
necessary  for  specifying  the  vibration  standards. 

Before  the  1970' s,  most  of  the  induction  of  the  ambient 
vibration  data  applied  the  "peak  value  inclusion  method".  On  a 
piece  of  graph  paper  with  value  of  vibration  versus  frequency, 
the  "peak  value  inclusion  method"  picks  a  curve  which  contains 
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most  of  the  data  points  while  rejecting  a  minority  of  isolated 
data  points  which  stuck  out  as  the  standard  vibration  curve'' 

This  method  is  very  easy  to  do  but  its  drawback  is  that  its 
inductive  result  is  too  crude  and  conservative.  At  the  same  time 
it  is  also  difficult  to  reasonably  specify  the  time  of  the 
vibration  experiment^- ^ . 

Based  on  a  series  of  research  into  the  various  data 
inductive  methods  available  to  date,  the  present  paper  proposes 
the  inductive  method  for  the  experimental  data  collected  from  the 
aircraft  ambient  vibration  experiment  -  the  state  inductive 
method.  The  results  obtained  by  using  this  inductive  method  can 
better  reflect  the  synthesized  vibration  characteristics,  and  it 
is  relatively  non-conservative.  In  addition,  this  method 
introduces  the  view  point  of  the  limited  life  span  design.  In 
other  words,  the  time  of  the  function  experiment  and  the  time  of 
the  extended  experiment  are  determined  according  to  the  actual 
time  duration  of  the  various  states.  This  makes  the 
specification  of  the  time  of  vibration  experiment  more  reasonable 
and  reliable.  This  is  why  the  state  inductive  method  is  more 
advanced  and  reasonable  when  compared  with  the  peak  value 
inclusion  method. 

It  is  necessary  for  us  to  know  the  state-time  spectrum  of 
the  aircraft  ambient  vibration  in  order  to  apply  the  state 
inductive  method.  This  is  to  say  that  we  have  to  know  the  times 
taken  up  by  the  various  typical  states  of  the  aircraft  during 
each  average  flight. 

The  state-time  spectrum  of  the  aircraft  ambient  vibration  is 
an  important  statistical  parameter.  Similar  to  the  loading 
fatigue  spectrum  of  the  aircraft,  it  is  obtained  statistically  by 
using  the  actual  flight  data  of  the  aircraft.  There  are, 
however,  major  differences  between  the  two:  1.  the  loading 
fatigue  spectrum  gives  the  repeated  number  of  every  degrees  of 
overloading  values  during  the  application  life  span  of  the 
aircraft,  while  the  state-time  spectrum  gives  the  times 


occupied  by  every  type  of  typical  states  (such  as  taxiing  during 
take-off,  climbing,  cruising,  diving,  and  spiraling)  during  each 
average  flight.  2.  The  compilation  of  the  loading  fatigue 
spectrum  does  not  take  the  conditions  of  the  ground  tests  into 
consideration.  However,  the  vibration  caused  by  the  ground  tests 
can  induce  the  same  kinds  of  damage  and  destruction  to  the 
products.  The  actual  measured  data  indicates  that  for  certain 
parts  of  the  aircraft  the  value  of  vibration  caused  by  the  ground 
tests  was  even  greater  than  that  produced  during  the  flight. 

This  is  why  the  construction  of  the  state-time  spectrum  of 
aircraft  ambient  vibration  will  have  to  simultaneously  consider 
both  the  conditions  in  the  air  as  well  as  the  ground  test 
conditions . 

The  aircraft  is  an  aviation  vehicle  that  will  be  used 
repeatedly.  There  are  many  different  kinds  of  flight  courses  and 
the  flight  states  are  also  complicated.  This  paper  presents  a 
method  for  the  compilation  of  the  state-time  spectrum  of  aircraft 
ambient  vibration.  This  method  is  established  on  the  foundation 
of  large  amounts  of  flight  statistical  data,  and -it  also  fully 
considers  both  the  flight  condition  in  the  air  as  well  as  the 
ground  test  condition.  The  stated-time  spectrum  of  aircraft 
ambient  vibration  compiled  this  way  can  better  reflect  the  actual  / 26 
application  condition  of  the  aircraft.  We  will  now  introduce  the 
compiling  method  as  follows. 

II.  COMPILING  PROCEDURES 

1.  After  checking  the  flight  records  of  the  aircrafts 
involved  in  the  statistics,  the  following  data  are  computed 
statistically: 

(1)  the  actual  numbers  of  take  of fs  and  landings  of 

each  aircraft  for  each  course  N..  and  the 

ij 

corresponding  flight  time  . 

(2)  the  numbers  of  take  of  fs  and  landings  of  all 
the  aircrafts  involved  in  the  statistics 
for  each  course  and  the  flight  time  Tj 


rTw^ 


N>-  S 

i  »  l 

2  7*;-, 

I  »  L 


(1) 

(2) 


(3)  the  actual  numbers  of  total  take  offs  and 
landing  of  all  the  aircrafts  involved  in 
the  statistics  N  and  the  total  flight  time 
T: 

N~  2  (3) 

i  •  i , 

r-  2  r,  (4) 

/  •  i 


2.  After  checking  the  record  of  the  propeller  of  the 
aircraft  (if  an  aircraft  has  two  or  more  propellers,  only  the 
record  of  one  of  the  propellers  will  be  included  in  the 
statistics),  the  ratio  between  the  flight  time  for  an  average 
flight  (including  take  off,  landing  and  taxiing)  and  the  time  for 
the  ground  test  is  statistically  computed  to  be  A:1. 

In  order  to  assure  statistical  accuracy,  the  numbers  of 
aircrafts  and  propellers  generally  should  not  be  less  than  30. 

3.  The  time  required  to  complete  an  average  flight  is 
computed  as 


4.  The  combined  classification  of  the  flight  courses: 
there  are  many  different  kinds  of  flight  courses  during  the 
actual  training.  However,  some  of  the  contents  of  these  courses 
are  the  same  or  similar  to  one  another  from  the  view  point  of 
ambient  vibration,so  they  can  be  combined.  The  principles  used 
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to  combine  these  courses  are:  the  corresponding  flight  time, 
flight  altitude,  flight  velocity  and  the  flight  motion  should  be 
roughly  similar  for  courses  which  are  combined  into  the  same 
category.  After  the  combined  classification,  the  course  with  the 
highest  numbers  of  actual  flights  in  each  category  will  be  used 
to  represent  a  typical  flight  course  in  the  same  category. 

By  adding  the  flight  times  of  the  courses  in  the  same 
category,  we  can  compute  statistically  the  flight  time  of  the  K- 
th  category  courses  for  all  the  aircrafts  involved  in  the 
statistics  T^. 

Assume  that  the  application  life  span  of  the  aircraft  is  T  , 

s 

the  flight  time  for  each  category  within  T  is 

s 

7'jx*“ Tt  (A.’**  1.  '-«••••  P)  (6) 

We  can  convert  T  ^  into  the  take  off  and  landing  numbers  for 
the  typical  flight  course  of  that  particular  category  as 

NtKmTs*/t, i  (A'»  1.  2 ,  •••,  p)  (7) 

5.  Two  methods  can  be  used  to  plot  out  the  state- time 
cross-section  of  a  typical  flight  course:  one  way  is  to  conduct 
test  flight  experiment  for  each  typical  course,  measure  the  time 
records  of  the  flight  parameters  such  as  the  flight  altitude, 
velocity,  flight  posture,  pressure  and  longitudinal  loading 
throughout  the  entire  flight  process  from  take  off  taxiing  to 
landing  taxiing;  the  other  method  is  to  conduct  theoretical 
estimation.  It  does  not  matter  whether  we  apply  one  method  over 
the  other,  or  if  we  use  both  methods,  the  objective  is  to  measure  / 2 7 
or  estimate  the  time  occupied  by  each  flight  state  in  each 
typical  course  as  well  as  the  corresponding  flight  altitude, 
velocity,  pressure,  and  the  longitudinal  loading  value.  In  this 
way,  we  can  obtain  the  state-time  cross-section  of  each  typical 


course 


6.  The  combination  of  the  flight  states:  the  aircraft  will 
have  to  perform  a  lot  of  flight  motions  during  the  actual  flight 
process.  At  the  same  time,  due  to  the  differences  in  the  state 
of  the  aircraft's  propeller  and  flight  parameters  such  as  flight 
altitude,  velocity,  and  longitudinal  loading  value,  the  flight 
state  of  the  aircraft  will  change  at  all  times.  It  is  necessary 
for  us  to  combine  and  classify  them  in  order  to  proceed  with  the 
statistical  analysis.  The  principles  for  the  classification  are: 
(1)  the  flight  motions  of  the  flight  states  in  the  same  category 
will  have  to  be  the  same  or  similar;  (2)  the  flight  altitude  will 
be  divided  into  low,  medium  and  high  altitude  ranges,  the  flight 
altitudes  of  the  flight  states  in  the  same  category  will  have  to 
be  within  the  same  altitude  range;  (3)  the  flight  velocity  will 
be  divided  into  3  to  5  levels  according  to  the  magnitude  of  the 
velocity,  the  flight  velocities  of  the  flight  states  in  the  same 
category  will  have  to  be  within  the  range  of  the  same  level. 

7.  By  adding  the  flight  times  of  the  various  states  which 
belong  to  the  same  state  category  and  the  same  course  category, 
we  can  obtain  the  time  occupied  by  each  category  -of  flight  state 
for  each  category  of  course  T^  (K  =  1,2,...,p,  Q  =  1,2,...r). 

3.  We  can  compute  the  time  occupied  by  each  flight  state 

category  within  T  as 

s 

7>  £  TkoNsk  (8) 

*  «  i 

9.  According  to  the  propeller's  ground  test  curve  we  can 
obtain  the  time  occupied  by  each  test  state  (such  as  slow  rate, 
specified  rate,  maximum  and  acceleration)  per  each  ground  test 
process  (a  *  1,2,...,  c). 
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Table  1.  An  example  of  the  state-time  spectrum  of  the 
ambient  vibration  of  the  aircraft 

2.  serial  number 

3.  air  space 

4.  low  altitude  air  space 

5.  medium  altitude  air  space 

6.  names  of  the  typical  states 

7.  take  off  taxiing 

8.  landing  taxiing 

9.  climbing 

10.  spiraling 

11.  diving  and  rapid  climbing 

12.  turning  and  increase  altitude 

13.  gliding  downward 


is 


14. 

15. 

16. 

17. 

18. 

19. 

20. 
21. 
22. 

23. 

24. 

25. 

26. 

27. 

28. 

29. 

30. 

31. 

32. 

33. 

34. 

35. 

36. 

37. 

38. 

39. 

40. 

41. 

42. 

43. 

44. 

45. 

46. 

47. 

48. 

49. 

50. 

51. 

52. 

53. 

54. 

We 


cruising  velocity  1 

cruising  velocity  2 

cruising  velocity  3 

put  down  the  landing  gear 

put  down  the  wing  flaps 

climbing 

spiraling 

diving  and  rapid  climbing 
turning  and  increase  altitude 
gliding  downward 
horizontal  rolling 
cruising  velocity  1 
cruising  velocity  2 

the  time  used  for  each  average  flight  (minutes) 
serial  number 
air  space 

medium  altitude  air  space 
high  altitude  air  space 
ground  tests 

names  of  the  typical  states 
cruising  velocity  3 
cruising  velocity  4 
climbing 
spiraling 

diving  and  rapid  climbing 
turning  and  increase  altitude 
gliding  downward 
horizontal  rolling 
cruising  velocity  1 
cruising  velocity  2 
slow  rate 

open  the  gas  release  belt 
close  the  gas  release  belt 
right  propeller  at  constant  rate 
right  propeller  at  maximum  rate 
right  propeller  thrust 
left  propeller  at  constant  rate 
left  propeller  at  maximum  rate 
left  propeller  thrust 

property  of  acceleration  or  deceleration 
the  time  used  per  each  average  flight  (minute) 

already  know  that  the  time  required  for  each  ground  test  / 28 
the  necessary  numbers  of  ground  tests  within  Tg  is  thus 

(9) 
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The  times  in  which  the  aircraft  is  in  the  various  test 

states  within  T  are 
s 

T.-V,a  c)  (10) 

10.  We  can  compute  the  average  numbers  of  take  off  and 
landing  within  T  as 

w 


(11) 


11.  The  time  occupied  by  each  state  for  each  average  flight 
process  is  thus 


(<?-  l,  2,-,  P) 
f.=»r./iV.  (o-i.  2.-.  c) 


(12) 


(13) 


This  is  the  required  state-time  spectrum  of  the  ambient 
vibration  of  the  aircraft. 

Table  1  shows  one  example  of  the  state-time  spectrum  of  the 
ambient  vibration  of  the  aircraft. 

III.  EXAMPLES  OF  APPLICATION 

1.  If  we  already  know  the  application  life  span  T  of  the 

s 

aircraft  or  the  on  board  equipment,  the  state-time  spectrum  of 

the  ambient  vibration  of  the  aircraft,  as  well  as  the  random 

2 

vibration  quantity  Wq  (Q  =  1,2,...,r)  (g  / Hz)  experimentally  for 
this  particular  equipment  for  each  state,  we  can  determine  the 
measured  value  and  time  of  the  vibration  experiment  corresponding 
to  that  particular  equipment. 

Reference  [13  specifies  that  the  maximum  W  among  all  the 

max  ° 

WQ  be  selected  as  the  measured  value  for  the  function  experiment. 
It  also  specifies  that  we  should  select  a  state  corresponding  to 
Wmax  ^rom  t^e  state-time  spectrum  and  set  the  time  of  an  actual 
flight  at  that  particular  state. 
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If  we  set  the  measured  value  for  the  extended  experiment  to 
be  WM,  as  well,  the  time  required  for  the  extended  experiment 
per  each  average  flight  is 


t.- 


(14) 


The  take  off  and  landing  numbers  of  the  aircraft  within  T 


are 


+- i-)/f  (15) 


The  time  of  the  extended  duration  experiment  necessary  for  that 
particular  equipment  is  thus 


2.  To  determine  the  vibrational  life  span  of  the  on  board 
equipment!  a  certain  piece  of  equipment  will  suffer  unacceptable 
damage  or  destruction  after  it  has  been  subjected  to  a  vibration 
quantity  W  for  a  time  of  T  .  Now  we  will  determine  the 
vibrational  life  span  of  that  particular  piece  of  equipment. 

Based  on  the  state-time  spectrum  of  the  ambient  vibration  of 
the  aircraft,  we  can  obtain  the  time  te  necessary  to  conduct  an 
extended  duration  experiment  for  that  particular  piece  of 
equipment  per  each  average  flight  from  equation  (14). 

In  this  way,  the  vibrational  life  span  of  that  piece  of  on 
board  equipment  can  be  determined  by  the  equation  below: 


7\- 


(17) 


As  we  all  know,  the  application  life  span  of  the  on  board 
equipment  is  determined  by  many  factors.  In  addition  to  the 
fatigue  damage  induced  by  vibration,  factors  such  as  wearing, 
corrosion  and  mould  can  all  cause  the  equipment  to  lose  the 
functional  capability.  If  the  vibrational  life  span  T^  computed 
from  equation  (17)  is  smaller  than  the  application  life  span 


computed  based  on  other  factors,  we  should  use  the  vibrational 
life  span  Ty  to  determine  the  replacement  period  of  that 
particular  piece  of  on  board  equipment. 
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.  COMPILATION  AND  APPLICATION  OF  A 
STATE-TIME  SPECTRUM  OF  AIRCRAFT 
AMBIENT  VIBRATION  • 

Gong  Qingxiang 

(Af anchang  Aircraft  Manu factoring  Company ) 

Abstract 

A  state-time  spectrum  of  aircraft  ambient  vibration  is  an  important 
statistical  parameter  of  an  aircraft.  It  is  an  essential  precondition  for  ii- 
ducing  ambient  vibration  data  of  an  aircraft'  with  the  state  inductive 
method. 

This  paper  puts  forward  a  compiling  method  of  the  state-time  ^spec- 
trum  of  aircraft  ambient  vibration.  It  is  based  on  -a  great  deal  of  statis¬ 
tical  data  taking  both  flight  and  ground  tests  into  comprehensive  consi¬ 
deration.  Therefore,  the  state-time  spectrum  of  aircraft  ambient  vibration 
drav.  n  up  by  this  method  can  give  true  expression  of  the  aircraft  practice. 

This  method  has  found  engineering  application. 
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